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Cosmological observables are used to construct cosmological models. Since cosmological obser- 
vations are limited to the light cone, a fixed number of observables (even measured to arbitrary 
accuracy) may not uniquely determine a cosmological model without additional assumptions or 
considerations. A prescription for constructing a spherically symmetric, inhomogeneous cosmologi- 
cal model that exactly reproduces the luminosity-distance as a function of redshift and the light-cone 
mass density as a function of redshift of a ACDM model is employed to gain insight into how an 
inhomogeneous cosmological model might mimic dark energy models. 



I. INTRODUCTION 



>- 1 PACS numbers: 98.70.Cq 

O 

In the last decade or so remarkable progress has been made in measuring cosmological parameters to unprecedented 
r \ ' accuracy. Parameters such as the Hubble constant, Hq — 72 ± 8 km s _1 Mpc -1 1], the temperature of the cosmic 
background radiation (CBR), T = 2.728 ± 0.004 K 0], and many other parameters, are now known to impressive 
■ precision. Of course, we don't invest so much time and effort determining cosmological parameters because of an 
interest in numerology, but rather, because the parameters are necessary inputs for the task of constructing a standard 
Q ■ cosmological model. 

In turn, we do not construct cosmological models simply to make predictions to compare with observations, but 
rather, to guide us to, in the words of Einstein 0], "a deeper and more consistent comprehension of the physical and 
astronomical facts." 

The precision cosmological measurements have lead to the latest cosmological model, usually called the standard 
cosmological model, the concordance model, or simply ACDM, where A indicates the inclusion of Einstein's cosmo- 
logical constant (or more generally, dark energy), and CDM stands for cold dark matter. The model is but one 
possible realization of homogeneous and isotropic cosmological solutions to Einstein's equations. (We will refer to 
homogeneous/isotropic models as Friedmann-Lemai'tre-Robertson-Walker (FLRW) models.) The question yet to 
be answered is whether today's standard model will lead us to a deeper and more consistent comprehension of the 
physical and astronomical facts. 

A troubling feature of the ACDM model is the present composition of the universe: 95% is dark! Of the dark 
components, roughly 25% of the total mass-energy density is dark matter, associated with galaxies, clusters of galaxies, 
and other bound structures. The bulk of the mass-energy of the universe, about 70% in the standard ACDM model, 
is attributed to dark energy, which is usually said to drive an accelerated expansion of the Universe. 

The observations and phenomena that lead to the assumptions of dark matter and dark energy are real, but that 
does not necessarily imply that dark matter and dark energy as usually envisioned are real. Dark matter and dark 
!— 3 ' energy are good names for the phenomena. But naming is not explaining! 

In some sense, an even more remarkable feature of the standard cosmological model is that it seems capable of 
accounting for all cosmological observations; i.e., it seems to work! Here, we address the issue of whether the agreement 
with observations prejudices our judgment regarding the likelihood that the ACDM model is the final answer. In 
particular, we will discuss the possibility that dark energy, per se, does not exist. 

The most important point regarding dark energy (and also, for that matter, the acceleration of the expansion of 
the universe), which is central to the motivation for this paper, is that all evidence for dark energy is indirect. The 
only effect of dark energy is on the expansion history of the Universe. 

The expansion rate of the Universe is perhaps the most fundamental quantity in cosmology. In the standard 
cosmological model the expansion rate is determined from the 00-component of the Einstein equations, Goo = k Tqq, 
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where k = 8irG. In FLRW models, Goo = 3H 2 + k/a 2 , where H = a/a is the expansion rate (here, a is the Robertson- 
Walker scale factor and k = ±1 or 0, depending on the geometry). With the assumption of a perfect-fluid stress 
tensor, Too = where p is the mass density. 

The stress tensor is assumed to consist of several components with different equations of state. For any component 
i, the equation of state for that component is defined as Wi — Pi/ pi- If Wi is constant, the energy density in component 
i evolves with redshift z = oo/a — 1 as pi(z) = p;(0) (1 + z) zl " 1+Wz \ where ao is the present value of the scale factor. 
Then the expansion rate as a function of redshift can be expressed as 



H 2 (z) = Hi { ft A exp 



' IZl -3[l + w A (*i)] 



l + Zl 



+ n k (i + z f + n M (i + + + *f >■ (i) 



Here, tti is the ratio of the present energy density in component i compared to the critical density pc — 3H 2 /n. The 
subscript "fc" indicates a spatial-curvature contribution (wk = —1/3), "M" indicates a matter component (% = 0), 
"i?" indicates a radiation component (wr — 1/3), and "A" represents dark energy. In the event that the function 
w\(z) in the last term is independent of z, the last term in H 2 (z) becomes + Z ) 3 ( 1+WA \ The equation of state 
for a cosmological constant is wa = — 1. 

It appears that the expansion rate as a function of z depends on four constants {fifc, 17m, Qr, ^a} and one function 
wa(z). In reality, the situation is much more predictive. First, since H 2 (0) — H 2 , it follows that J^- fit = 1, removing 
one constant. 

The value of fl k is well determined by WMAP to be small, Q k = -0.0026tg;$^ at the 68% confidence level Q, 
effectively removing yet another constant. The value of _r is also well determined by CBR measurements (5 x 10 -5 ) 
removing a further constant. One may then make the less secure assumption that f2jv/ is "known," or more profitably, 
use observational constraints to marginalize over the fit's, leaving only the function wa(z) undetermined. The function 
must be determined through its effect on the expansion history of the Universe, which is manifest through various 
cosmological observables. 

Many cosmological observables depend on the expansion history of the Universe. They often depend on the expan- 
sion history through the coordinate distance r of a source of redshift z. Recall that the Robertson- Walker metric can 
be written in the form 



ds 2 = dt 2 ~a 2 (t) 



dr 2 



1 — kr 2 



r 2 dVL 2 



(2) 



where here d£l is the angular differential and r is the comoving "radial" coordinate. The radial coordinate of a source 
at redshift z is determined by integrating the null geodesic equation (ds 2 = 0) to obtain 




dz 1 
H(z') 



(3) 



where sin, 1, sinh obtains for k — +1, 0, 1, respectively. 

Observables such as the luminosity distance, the angular-diameter distance, the volume element, and the age of the 
Universe, all depend on the time evolution of the expansion rate, hence on the properties of the stress tensor, hence 
on the dark energy fluid characterized by wa ■ 

At present, a constant wa{z) seems to fit the data, and furthermore wa = —1, the value of wa for a cosmological 
constant, is not disallowed. It is also true that since there is no apparent reason to exclude vacuum energy from the 
stress tensor, there is no reason not to imagine it enters the determination of the expansion rate. The problem is 
that the magnitude of the expected vacuum energy is very many orders of magnitude larger than expected. If the 
cosmological constant is part of nature, then either some principle must be found to explain the much smaller-than- 
expected magnitude of 17a, or one must resort to anthropic arguments. 

This situation has led many to consider alternate explanations of the observations: modified gravity, new long-range 
forces, new ultra-light scalar fields, extra dimensions, branes, bulk, Lorentz-invariance violation, etc. Some of these 
approaches are, in many respects, more drastic than a small cosmological constant. For instance, theories of modified 
gravity typically do not respect the principle of general covariance. Since we should not abandon our principles 
without good cause, it is tempting to retreat to the position that since a cosmological constant is adequate to account 
for the data, just assume that it exists and call it a day? 

We resist the retreat because our goal in cosmology is more than just explaining the observations. The ingredients of 
a cosmological model must be deeply grounded in fundamental physics. Dark matter, dark energy, modified gravity, 
mysterious new forces and particles, etc., unless part of an overarching model of nature, should not be part of a 
cosmological model. We may plant new ideas and features into cosmological models, but they will ultimately prove 
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barren unless grounded in fundamental physics. No matter how well a concordance model fits the data, its ingredients 
must not be discordant with fundamental physics. 

In this paper we explore a less traveled path, and consider the possibility that the Friedmann equation is not an 
adequate description of our inhomogeneous Universe [f| U 0, H, M, OH [HI E2) tHi EH EH • In this approach there 
is no need for dark energy, so no need for new long-range forces, modifications of general relativity, new ultralight 
particles, or anthropic reasoning. The phenomenon of inhomogeneities mimicking dark energy is sometimes referred 
to as "backreaction." 

In this paper we follow the prescription of Mustapha, Hellaby, and Ellis [l6| [MHE], recently applied to ACDM 
models by Celerier, Bolejko, Krasinski, and Hellaby [17], 1 and explore the implications of the fact that one can 
construct a spherically symmetric inhomogeneous model that exactly reproduces the redshift dependence of the 
luminosity distance and mass density of any given ACDM model. We show that in general one expects the cosmological 
solution thus obtained to have a "mixmaster" approach to the initial singularity, with the radial scale factor diverging 
and the angular scale factor vanishing. We then show that there is another cosmological solution with the same 
density and curvature functions as the singular model but with regular initial conditions where both scale factors 
vanish at the singularity, and that this regular solution well approximates another ACDM model. We then discuss 
the utility of different averaging prescriptions for understanding the results. 

In the next section we briefly review spherically symmetric cosmological solutions. In Sect. IIIII we review the 
construction of such an inhomogeneous model that reproduces key observational features of ACDM. In Sect. IIVI we 
evolve the model off the light cone, paying special attention to its behavior near the bang. In Sect. [V] we calculate 
some spatially averaged features of the model. Section IVT1 presents the conclusions and discusses their implications. 



II. LEMAITRE— TOLMAN— BONDI MODELS 



In this section we discuss Lemaitre-Tolman-Bondi 120, 2JJ, [22| (LTB) models in general, and in the next section we 
discuss a particular LTB model that exactly reproduces select observational features of the ACDM model. 

LTB models are spherically symmetric cosmological solutions to the Einstein equations with a dust stress-energy 
tensor. The LTB model is based on the assumptions that the system has purely radial motion and the motion is 
geodesic without shell crossing (otherwise the pressure could not be neglected). The line element in the synchronous 
gauge can be written in the form 

ds 2 = -dt 2 + R ' 2 ^\ dr 2 + R 2 (r, t) dil 2 . (4) 
1 + 8(r) 

Here, the prime superscript denotes d/dr. We will denote d/dt by an overdot. The function 8(r) is an arbitrary 
function of r. The Robertson- Walker metric can be recovered if R(r, t) — > a(t)r and f3(r) — ■> —kr 2 . 

In spherically symmetric models, in general there are two expansion rates: an angular expansion rate, H± = 
R(r,t)/R(r,t), and a radial expansion rate, H r = R'(r,t)/R'(r,t). With the dust equation of state, the Einstein 
equations may be expressed as 

8(r) B'(r) 
Hl(r,t) + 2HAr,t)HJr. t) - J*L. - = sp(r , () 

6 wt) + 2g!(r '" ~ 2 mk - m ' {r - t)H ^ " + wx) = -"*'•«• (5) 

These represent the generalization of the Friedmann equation for a homogeneous/isotropic universe to a spherically 
symmetric inhomogeneous universe. 

The equations may be manipulated to result in an easily integrable dynamical equation for R(r,t): 



Similar formalism was developed by Chung and Romano 18] to fit only the luminosity-distance, and by Yoo, Kai, and Nakao Il9t who 
simultaneous fit the luminosity distance and require no spatial variations of the age of the Universe. 
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where a(r) is a function of r related to the density p(r,t) by 



^ (r '* )= ^M)2(r,t)- (7) 
Equation © can be differentiated to yield the dynamical equation for R'(r,t): 

ft, t) = 0V) + a'(r)/R(r 7 1) - a{r)R'{r, t)/R 2 (r, t) 

1 ' ' 2R(r, t) ' K ' 

The two functions a(r) and [3(r) define the LTB model, playing roles analogous to the total energy density parameter 
f2 and the spatial curvature parameter f2fc of FLRW models. To see this, recall that the evolution of the scale factor 
and time for a hyperbolic dust Friedmann model evolves as (77 is conformal time) 



a(t) = a Q — — [cosh 77 - 1] , t-t BB =H — — [sinh r\ - 77] , (9) 
2£2fc 20,7 



where ao is the present value of the scale factor and tBB is the time of the big bang (usually set to zero). 

Since Friedmann models are homogeneous and isotropic, they are zero-dimensional models. 2 The LTB models are 
spherically symmetric, so they are one-dimensional problems. So for LTB models the density parameter Slo is replaced 
by a density function a(r) and the curvature parameter is replaced by a curvature function (3(r). If (3{r) > 0, 
which is the analog of hyperbolic Friedmann models, the evolution of the metric function R(r, t) can be described by 

R{r, t) = [cosh ry(r) - 1] ; t - t BB (r) = [sinh ^(r) - V (r)} , (10) 

where now, in general, the bang time, tBB{r) is a function of r, as is the effective conformal time rj{r). 

The photon geodesic equation for the position of the photon as a function of time, t(r), is found from the LTB 
metric: 



The redshift of the photon, z(r), is 



di(r) _ R'(r,i(r)) 



Equations (TTT]) and (JT2J) are solved with initial conditions z(r = 0) = and t(r = 0) = 0. The luminosity distance is 
then simply 

d L (z) = (l + z) 2 R(r,t(r)). (13) 
(We denote the luminosity distance by d^ to emphasize that it is a light-cone observable.) 

III. RECONSTRUCTING AN LTB MODEL WITH ACDM OBSERVATIONAL FEATURES 

In this section we will review the procedure and results that allow one to construct an LTB model that reproduces 

1. the luminosity-distance-redshift relationship, dh{z) 1 of a fiducial ACDM model, and 

2. the light-cone matter density as a function of redshift, p(z), of the fiducial ACDM model. 



2 The dimension of a model refers to the number of dynamically independent combinations of spatial dimensions that enter in physical 
quantities. 
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The "fiducial" ACDM model may in principle be any FLRW model with a cosmological constant. In this paper 
we will assume for the fiducial model a spatially- flat Universe, so Qm + = 1- Unless otherwise specified, we will 
also assume f^A = 0.7 (so VLm = 0.3). The assumption of spatial flatness will simplify our calculations. Neither that 
assumption, nor the exact value of Q\, is important in our considerations. 

The reconstruction method we follow was discussed in general by Mustapha, Hellaby, and Ellis in 1998 [l6| [MHE], 
and recently applied to ACDM models by Celerier, Bolejko, Krasinski, and Hellaby [CBKH], [l7j. Related formalism 
was developed by Chung and Romano [l8[ and by Yoo, Kai, and Nakao [l9j . 

Let us review the procedure. To follow the discussion, it is useful to refer to Fig. [1] Observers obtain information 
about the universe by means of photons, so our direct information about the Universe is limited to the light cone 
given by Eq. (fTTj) . Quantities evaluated on the light cone are denoted by a "hat." Therefore, the metric functions 
R(r,t) and R'(r,t), and the density p(r,t) evaluated on the light cone, are 

R(r,t(r)) = R; R'(r,t(r)) = R'; p(r,t(r))=p. (14) 

The light cone for an observer situated at the origin is illustrated in Fig. Q] by the solid line. The shaded area 
indicates the causal volume of the Universe within the light cone. Again, we emphasize that all direct astronomical 
information is limited to the light cone. We often speak of cosmology within our causal volume, or perhaps the density 
field at some fixed time like the present (which includes an acausal region of the Universe). We do this even though 
we have no direct observational evidence for it. To speak of the Universe in this way we must employ a cosmological 
model that allows us to depart from the light cone (in theory, not literally). What we say about the Universe off 
the light cone is only as trustworthy as the cosmological model used to generate the information. But as emphasized 
by MHE, light-cone observations (or at least the two light-cone observations we use) do not uniquely determine the 
cosmological model on the light-cone (i.e., both the ACDM and LTB models give the same results), so of course they 
do not uniquely determine the cosmological model off of the light cone. 

As in MHE, we take advantage of a coordinate freedom to simplify the calculation. The radial coordinate r has no 
physical significance; it can be rescaled in any convenient way. Here, we rescale r such that on the light cone 

R' = HoWl + P(r). (15) 

In order to find a(r) and f3(r) we follow the procedure outlined by MHE and followed in CBKH. Here we take the 
view that we will construct an LTB model (i.e., find a(r) and (3(r)) that agrees with observations of di,{z) and p(z) 
only out to r corresponding to z — 2, and beyond that we have no useful information about oIl(z) or p(z) that would 
constrain a(r) and (3(r). Perhaps information from cosmic microwave background (CMB) measurements could fill 
that niche, but since the purpose of this paper is to investigate constraints from the local Universe, we will not pursue 
constructing the model beyond a value of r corresponding to z = 2. 

1. The first step is to use the assumption that the LTB luminosity-distance-redshift relationship matches that of 
the fiducial ACDM model. In spatially flat ACDM, the luminosity distance is given by 

d L (z)=(i+z) r dz \ (i6) 

Jo -Wacdm \Z\ ) 

where Hacbm(z) = -ffo\/^A/(l + z ) 3 + ^A- Therefore, using Eq. (fTB"]) . 

R(z) = f - ^ (17) 

2. The next step is enforce the assumption that the observed mass density on the light cone in the LTB model as 
a function of z, p(z), matches the observed mass density as a function of z in the fiducial ACDM model. Using 
p(z) d 3 VLTB = PM.AGDM d 3 V AC DM, we find 3 



2 



H^R 2 {z)Kp{z) — = Sl M Z 



dz\ 



dz Hacdm(z) 
where we have used the fact that Pm,acdm(z) — QmSHqH^ 1 ^ + z) 



Hact>m{zi) 



(18) 



3 In MHE and CBKH, H~ 1 R 2 (z)Kp LTB (z)dr = nfhndz. 
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FIG. 1: This figure illustrates features and notation of the LTB model analyzed in this paper. The time coordinate is defined 
such that its value at the present time is to; to is not the present age of the Universe. In LTB models the age of the universe 
is a function of r, since, in general, the bang time, tsB, is a function of r. The age as a function of r is given by to — tss(r). 
The light cone (the past-directed null geodesic) for an observer at r = is denoted by the solid line. A "hat" indicates that a 
quantity is evaluated on the light cone, and hence is a function of only one variable: f(r) = f(r,t(r)). Hatted quantities are 
often written as functions of z, since that is the light-cone coordinate accessible to observations. The r coordinate is scaled 
such that t(r) is given by t(r) — to — H^ l r. The values of R and R' on the light cone are denoted by R = R(r,t(r)) and 
R' = R'(r,t(r)). The shaded area is our causal past. 



3. We make use of the fact that 
H dR(z) 



dz 
dr 



dz 



(1 + *) 



1 f z — dr 

- / Ho 1 Kp(z)R(z 1 )(l + z 1 )—dz 1 
l J az\ 



= (l + z) 



^acdm(^) 
H • 



(19) 



where the last equality (not found in MHE or CBKH) obtains only for a spatially flat Friedmann model. The 
resulting z(r) is shown in Fig. [2] 



4. We then solve the differential equation from MHE for a(r): 

da 1 i c-o , , > I 

— = -H Q 1 R 2 (z) K p(z) 



HodR/dr 



1 - 



R 



Ho 



dR 

dr 



(20) 



where dR/dr = dR/dz ■ dz/dr. A The initial condition is cv(0) = 0. 
5. Finally, /3(r) is given by 

/ v 2 

da 



1 



dr H ~ 1 R 2 Kp / 



(21) 



The values of a, (3, a', and /?' obtained using this procedure are shown in Fig. [3] They agree with the results 
of CBKH. 



4 Here we mention a subtlety in Eq. il 2 ( ) t) . The factor dR/dr vanishes at some value of z (in the fiducial ACDM model, at z ~ 1.6), 
but since 1 — a/R also vanishes at that point, the equation is regular. However, care must be taken when numerically integrating the 
equation. 
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FIG. 2: The redshift z and the present proper distance Hodp(to) as a function of the coordinate r. 
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FIG. 3: The functions a(z), /3(z), ot (z), and /3'(z) as a function of z. 



We emphasize that for these values of a(r) and f3(r), the LTB model exactly reproduces the luminosity-distance- 
redshift relation and the mass-density-redshift relation of ACDM. 

Once a(r) and j3(r) are determined, along with R and R' , one may use Eq. ((8]) to solve for R'(r,t) and either Eq. 
© or Eq. (fit)]) to solve for R(r,t) from the bang time, iss(r), until today, to. 

Again, we return to a point made earlier: We only have observational information on our light cone. This means 
that although one might speak of the density as a function of r at a particular time, it is not accessible to observations. 
Nevertheless, it is instructive to examine p(r,t) at fixed times. Such an example is the present value of p, p(r, to). 
CBKH present it as a function of the metric function to), but here, we present it in terms of a physical quantity, 
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FIG. 4: The present value of p, p(to), as a function of the present proper distance. It is normalized to the value in the fiducial 
ACDM model, Pm-acdm 



namely dp(to), the present proper distance to r. The present proper distance is given by 



and is shown in Fig. [2] as a function of coordinate r. 

Again, as noted by CBKH, while previous attempts to construct an LTB model to reproduce dj,(z) of the fiducial 
ACDM model have resorted to a large local void % [H [H, [H, S3, El, S3, S3, M, IU M, H, EE &L Ezl, El , here we 
find a very mild increase (less than 4%) in p(r). (As already remarked, this increase is unobservable since only p(z) 
is observable, and that exactly agrees with ACDM). 

A point that will be developed further is worth introducing here. It has often been remarked that it is not possible 
for inhomogeneous cosmologies to reproduce cIl(z) of ACDM models because large, unobserved, inhomogeneities 
would be necessary. This is usually stated as either 1) The large inhomogeneities would result in large, unobserved, 
peculiar velocities, or 2) One can always write the metric in the perturbed Newtonian form. These two arguments 
are often said to exclude the possibility that an LTB model can mimic dark energy observables. 

However, the present consideration seems to provide (yet another) counterexample to that statement, since in the 
LTB model constructed here, cIl(z) is exactly that of the fiducial ACDM model while p(z) is also exactly that of the 
fiducial ACDM model. 

In the next section we explore the time evolution of the model off of the light cone. 



Let us first examine the evolution of R(r, t) and R'(r, t) for the model constructed in the previous section. Again, 
we only consider the evolution out to r corresponding to z = 2. The values of R(r, t) and R'(r, t) for several values of 
r (labeled by the corresponding values of z) are indicated in Fig. [5l 

The striking feature is that for r ^ 0, as tss(r) is approached, R'(r,t) — > oo. 5 (The solution at the origin at 



5 To better appreciate why this feature is striking, recall that for an FLRW solution, R'(r,t) — > a(t), where a is the scale factor. One 
would expect the scale factor to approach zero at the bang time. Here, the effective radial scale factor R! blows up as the angular scale 
factor R vanishes. 




(22) 



IV. FROM THE BANG TIME UNTIL TODAY 



9 



ACDM 




FIG. 5: The evolution of R(r,t) and R'(r,t) with time for the model discussed in Sec. Mil Labels on the curve correspond to 
values of r corresponding to the indicated values of z. R(0,t) — is not indicated. The age of the universe in the fiducial 
ACDM model is shown by the dotted vertical line. 



the bang time is R'(0, t BB (r)) = 0.) The behavior of the radial scale factor increasing as the angular scale factor 
decreasing is reminiscent of a type of Type I Bianchi model (the so-called Kasner solution) . 

To understand this behavior, consider the evolution of R' at early times. It is most convenient to work with the 
parameter r\ rather than time. From Eq. (jTDJ) , dt/d/q = R(r, i)//? 1 / 2 , and using Eq. (JSJ), we find 

dR' _ a' R' 
dr\ Aj3 rj 

Two small-77 solutions are possible. First, if one assumes that R' — > as rj — > 0, then the small-77 solution is 

R' -> — 7y 2 , (24) 
12/3 ' ' V ' 

which vanishes at the bang time (77 = 0). We will refer to this as the regular solution. 

The other solution to Eq. (j23|) obtains if as r\ — > 0, the second term on the rhs of the equation dominates. In this 
case, the solution is 

R! - -, (25) 
V 

where A is an constant. We will refer to this as the singular solution. 

The different behaviors near the origin can be understood by expressing R'(r,t) as [39j 

R'(r,t) = {^-j) R(r,t) - - [t - t BB {r)] R(r,t) - t' BB (r)R(r,t). (26) 

At* = t BB (r), R{r,t BB {r)) = 0, and Eq. {26]) implies R'(r, t BB {r)) = -t' BB (r)R(r,t BB (r)). Since R(r, t BB (r)) = +00 
and in our case t' BB (r) < 0, we see that R'(r, t BB (r)) = +00 as seen in Fig. [5j Since as t — > 0, 

M a(r) T] 3 (r)^ R ^^l^_. £ , 2 P 1/2 ( 2 7) 



/?3/2( r ) 12 ' (3{r) 4 ' 

the constant A in Eq. (gSJ) is identified to be A = -i^ B (r)2/3 1 / 2 
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FIG. 6: The evolution of R(r,t) and R'(r,t) with time for a(r) and /3(r) from the model discussed in Sec. IIIII but with the 
condition that t' BB (r) = 0. Labels on the curve correspond to values of r corresponding to the indicated values of z. 



The solution R'(r, fss(r)) = only will occur if t' BB (r) — 0, or in other words, if the big-bang surface is synchronous. 
In this case we can set tBB = 0, and one can easy recover Eq. (fM| from Eq. (j2l)|) . 

Different initial conditions, either Eq. (fM| or (j2"5|) , will lead to different evolutions of R'ijf) for the same values of 
R{rj). It is not difficult to see that the space spanned by singular solutions is larger than the space spanned by the 
regular solution because the constant A (related to t' BB (r)) of the singular solutions is undetermined. The regular 
initial condition leads to a unique value of R' at some later time, whereas the singular initial condition can lead to any 
value of R' at some later time, modulo the value obtained from the regular initial condition. Or as we have shown, 
the regular initial condition only occurs for a simultaneous bang surface. 

There is no particular reason to believe that R' would correspond to the unique value obtained from the regular 
initial condition. This conforms to the expectation that unless there is some symmetry to prevent divergence, the 
generic approach to the singularity should not be regular. 

At this point it should be pointed out that the small-77 behavior might be irrelevant for three reasons: 1) The LTB 
model is merely a toy model, not to be taken seriously. 2) The LTB model we consider has only a dust component 
to the stress tensor. Adding radiation may change the outcome, or imagining an early inflationary phase might also 
change the outcome. 3) Finally, perhaps it is premature to reject the singular solution without fully exploring why it 
is inappropriate, i.e., what is the observational basis for excluding it? 

We eschew these issues and ask a modest question of how the observations that went into the construction of a(r) 
and /3(r) would be changed if we substituted the regular initial conditions. This is quite straightforward: we have a(r) 
and /3(r), all we have to do is choose initial conditions R(r,0) = R'(r,0) = 0, evolve the metric functions according 
to Eqs. © and ©, enforce t' BB {r) — 0, and calculate the observables <Ll{z) and p(z). 

The results for R and R 1 following this procedure are illustrated in Fig. [5J 

Now consider the observational difference between the regular and singular solutions. In Fig.[7]we show the difference 
in di(jzt) for the two solutions. The regular solution is distinguishable from the singular solution, but the singular 
solution closely resembles another ACDM model, one with £7a = 0.5 rather than Oa = 0.7 of the fiducial model, at 
least out to redshift two. In fact, out to redshift unity, the difference is less than about 0.02 magnitudes. 

Just as the luminosity distance for the regular cosmological solution differs from the fiducial ACDM model, the 
density does as well. This is illustrated in Fig. [5] Note that the reference ACDM model for this figure is the model 
with 11a = 0.5, which well fits the dt,{z) for the regular LTB model. This figure well illustrates that the density along 
the light cone compared to a fiducial model can decrease as a function of z, while at the present time, it can increase 
as a function of r. 

It is beyond the purpose of this paper to see whether one can fiddle with a(r) and f3(r) to find a regular solution 
that is a better fit to the fiducial ACDM model. 
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FIG. 7: The luminosity-distance vs. redshift for several cosmological models. The regular LTB model is the LTB model with 
a(r) and (3(r) from Sect. Mil but with regular initial conditions. Also shown for comparison is <1l{z) for a ACDM model with 
S7a = 0.5 instead of the fiducial model of this paper (Oa = 0.7). The singular LTB solution is constructed to give the exact 
(Il(z) as for the fiducial ACDM model. Finally, in keeping with tradition, the result for the Einstein-de Sitter model (spatially 
flat dust cosmology) is shown. Here A(m — M) is the difference in magnitudes between the indicated model and the open, 
empty cosmological model. (M is a nuisance parameter related to the absolute magnitude of supernovae.) 
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FIG. 8: The density as a function of redshift on the light cone (solid curve) and at the present time (dashed curve) for the 
regular LTB model constructed in Sec. Mil For the regular LTB model, the reference ACDM model has = 0.5. 



The conclusion of this section is that if one attempts to construct an LTB model that agrees exactly with di,{z) 
and p(z) of the fiducial ACDM model, the resulting model will have singular initial conditions for R' . We have argued 
that this is to be expected. We surmise that this will be true in general, and is not specific to our choice of the fiducial 
ACDM model. We also demonstrated that one can find a regular solution with the same values of a(r) and (3(r) that 
results in a cIl{z) and p(z) that approximates a ACDM model close to the fiducial model at late time. 

In the next section we return to the singular models and examine whether they can be described in terms of averaged 
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quantities. 

V. DESCRIPTION IN TERMS OF AN AVERAGING PROCEDURE 

One proposed route to understanding the role of inhomogcncities in mimicking dark energy is by employing an 
averaging procedure. Our calculations allow us to study this proposal. First, let us review the averaging procedure. 
A fundamental quantity in the analysis is the velocity gradient tensor, which is defined as 

©'i = u\j = ^h ik h kj , (28) 

where h n is the spatial metric and is the fluid four velocity. The tensor Z • represents the extrinsic curvature of 
the spatial hypersurfaces orthogonal to the fluid flow. It may be decomposed in terms of a trace term and a traceless 
tensor as 

<->', <-M', • o' r (29) 

where O is the volume-expansion scalar. The traceless tensor a 1 ^ is the shear. 

The evolution equations for the expansion and the shear come from the space-space components of Einstein's 
equations (see e.g., Ref. Combining the expansion evolution equation with the energy constraint gives the 

Raychaudhuri equation, 

e + ie 2 + 2o- 2 + ~«p = o. (30) 

From the Raychaudhuri equation it is straightforward to verify that local fluid elements cannot undergo accelerated 
expansion. (This point was emphasized by Hirata and Seljak [4lj.) But, of course, as now appreciated, that point is 
irrelevant. Locally the expansion does not accelerate, but it is incorrect to assume that acceleration does not occur 
when the fluid is coarse-grained over a finite domain. The reason is trivial: the time derivative of the average of O 
and the average of the time derivative of O are not the same because of the time dependence of the coarse-graining 
volume. 

There has been a lot of work recently regarding the averaging procedure. First, we take the original proposal and 
average at a fixed time. We will then average over the light cone. 

Let us denote the coarse-grained value of a quantity T by its average over a spatial domain D: 



_J D VhJ 7 d 3 x ^ 

Following the work of Buchert [H, |43[ , we define a dimensionless scale factor 

a D (t) = (J^\ ' ; V D = j Vhd 3 x, (32) 

where Vn is the volume of our coarse-graining domain (the subscript "0" denotes the present time) . In the LTB model 
Vh = R 2 (r, t )R'(r, t )/y/l + /3(r). 

The coarse-grained Hubble rate Hp will be 

H D = ?° = \(Q) D . (33) 
an 3 



The smoothing procedure leads to the evolution equations for the coarse-grained scale factor [42J, |4 

— = -7 (Pcff + 3p e ff) , (34) 
a D 6 
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where the effective energy density and pressure terms are ( 3 TZ is the spatial curvature) 

Peft 
PeS 

and we have introduced the kinematical backreaction 

Q D = \{{Q 2 ) D -{Q? D )-2{a 2 ) D . (38) 

Finally, we can define an effective equation of state, u> e ff, to be w e fi = p e s/ p e s- 
Now for LTB models, 

r fe -2^ + ^ 

2E( r Sfe) -g(E r ofcJ =s{r^r:) 

0' 

~ 2 i~ 2 RK- < 39 > 

The various terms that enter the calculation of p e g and p e g are shown in Fig. [9] as a function of the proper distance 
that sets the scale of the averaging volume. Note that the shear is small, as is the kinematical back reaction since 
(Q 2 )d — Shown in Fig. [TO] is the effective value of w, w e g = p c s/PcS- Note that if one models the fiducial 

ACDM model as a single fluid and calculates w, the value today is w — —0.7. 

Another approach is light-cone averaging. Here, what is of interest isn't the average on a (acausal) constant-time 
hypersurface, but rather the average on a hypersurface coinciding with the light cone. One way to approach this 
is just to substitute the "hatted" quantities in Eqs. (|32j|39[) . The resulting terms that enter the calculation of the 
averaged quantities are shown in Fig. 1111 and the effective value of the equation of state parameter is shown in Fig. 
IT21 Also shown in Fig. [12] is the equation of state parameter for the fiducial ACDM model if one would model the 
total energy (including matter and A) as a single fluid. 

There are two things we learn from the exercise of this section. First, neither averaging on a constant-time hyper- 
surface nor light-cone averaging is easy to connect with the observations corresponding to parameters of the ACDM 
model. However, for both averaging prescriptions the effective equation-of-state parameter is negative, indicating 
acceleration of some averaged scale factor. This may be a valuable check in more complicated models. It is interest- 
ing to note that the shear contribution to the evolution of the averaged quantities is insignificant in both averaging 
prescriptions. This seems to be a counterexample to naive claims that if one constructs an inhomogeneous universe 
to mimic the observations in ACDM one would generate unacceptably large shear terms. 



(P) 



D - 

Qd 
2k 



Qd_ 
2k 

6k 



2k 



(36) 
(37) 



e = 

3 K = 



VI. DISCUSSION AND CONCLUSIONS 



In the discussion of what we have learned in this investigation, it is useful again to refer the reader to Fig. [JJ and to 
reiterate that all direct cosmological information we have is limited to the light cone, illustrated by the solid line in 
the figure. Any statement about the Universe off of the light cone is unsupported by direct observations, and hence 
requires a cosmological model for the evolution of the Universe. 

An interesting lesson from the consideration of this paper is that if one is just restricted to the two cosmological 
observables employed in this paper, di,(z) and p(z), it is impossible to determine uniquely a cosmological model 
fl6L In our case, the fiducial ACDM model and the singular LTB model result in identical observations for cIl(z) 
and p(z). That result is valid even if one imagines perfect astronomical observations, since by construction the two 
cosmological observables are degenerate in the two models. 

There are two ways to break the degeneracy, one observational and the other theoretical. One promising avenue 
to break the degeneracy with observational data might be to take advantage of the fact that in LTB models there 
are two different expansion rates, a radial expansion rate H r , and an angular expansion rate H± (they were defined 
in Sect. [TT|. Figure [TBI shows the differences between the expansion rates in the singular model. Note that the fact 
that H r = H\cdm is guaranteed by the input requirements used to construct the model. Quartin and Amendola [44j 
have recently discussed using cosmic parallax and redshift drift to distinguish between void models and dark energy. 
It is possible to pile on additional observational constraints to break the degeneracy. 
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FIG. 9: The various terms that enter the calculation for p c g and p e ff for the averaging at the present time as a function of 
proper distance that sets the scale of the averaging volume. 



I ' 1 ' 1 1 r 




FIG. 10: The value of w e s as a function of proper distance that sets the scale of the averaging volume for the averaging at the 
present time. 



While this can be done, it misses the important point that if one is restricted to light-cone observations, it is 
unclear whether any number of observations alone, even of arbitrary precision and with no systematic uncertainties, 
can uniquely determine a cosmological model, even if it is possible to rule out any LTB model. An FLRW model is 
a zero-dimensional (dimensions of space) model, and the LTB solution is a one-dimensional model. Intuition from 
other fields lead us to believe that the range of solutions in two-dimensional, or fully inhomogeneous three-dimensional 
models, will be richer. So it is unclear whether any number of cosmological observations, even to arbitrary accuracy, 
can uniquely specify a cosmological model. 

This point was discussed by Kolb, Marra, and Matarrese (KMM) (45[. In that work they introduce the concept 
of cosmological background solutions. This concept will be useful in elucidating several points, so we briefly review 
some concepts from KMM. 
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FIG. 11: The various terms that enter the calculation for p e ff and p e ff for light-cone averaging as a function of z. 




FIG. 12: The value of w e ff as a function of z for light-cone averaging. 



A cosmological background solution is denned as a mean-field geometry of suitably averaged Einstein equations, in 
which the average expansion is described by a single scale factor. A cosmological background solution depends upon 
the spatial curvature and the mass-energy content. Associated with the mass-energy content is a stress-energy tensor 
that describes a fluid (or fluids) with equation(s) of state that satisfy local energy conditions. 

KMM then define several types of background solutions; of relevance here is the Phenomenological Background 
Solution (PBS): a background solution that effectively describes the observations, that is, the data on the past light 
cone of an observer. The energy content and curvature of the PBS are not necessarily related to the spatial averages 
of the energy content and curvature of the observable universe. The equation of state of the stress-energy tensor of 
the PBS need not satisfy any of the local energy conditions of the mass-energy content of the universe. 

In this nomenclature, one might imagine that the true background solution is the LTB solution, or more generally 
some inhomogeneous model, and the fiducial ACDM solution is merely a phenomenological background solution. It 
is unclear whether it is possible, in principle, to exclude this possibility just on the basis of light-cone observations. 



16 



1.2 



1.1 



1 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

z 

FIG. 13: The expansion rates as a function of redshift compared to the result for the fiducial ACDM model. The expansion 
rates H r and H± were defined in Sect. [TT1 



Observations alone may not be able to determine the true cosmological model, but as remarked in the introduction, 
we expect the cosmological model to be part of our deeper comprehension of the physical and astronomical facts; it 
must fit into the larger framework of particle physics, general relativity, astrophysics, and cosmology. 

As an example of larger considerations, consider the theory of structure formation. In the fiducial ACDM model, 
structure grows in a background that is close to homogeneous and isotropic, but with small seed perturbations initially 
generated by inflation. Even in this zero-dimensional model, structure formation is a complicated affair (at least in 
the nonlinear regime) . Structure formation in LTB models is an even more complicated affair (23l . |24| : each spherical 
shell evolves as an independent FLRW model. One would expect the evolution of structure formation in LTB models 
to differ from the fiducial ACDM model. However, since structure formation in a three-dimensional model is more 
complicated still, it is difficult to see how structure formation might be conclusive. The same argument might apply 
to other cosmological tests like the integrated Sachs- Wolfe effect. 

Now, we return to the argument that purports to prove that inhomogeneous models cannot mimic the effects of 
dark energy (see, e.g., Refs. 0, |4?], EH]). The argument might be paraphrased as follows: Perform an expansion 
of the FLRW metric to Newtonian order in potential and velocity but take into account the fully nonlinear density 
inhomogeneities. In the conformal Newtonian gauge, the metric will be of the form 

ds 2 = a 2 (?7)[-(l + 2ip)dif + (1 - 2(j>)dx 2 }, (40) 

where rj is conformal time and <f) ~ ip. The volume average of the equation then yields [48| 

(jj) 2 = ~(p)(l-5W + 2K), (41) 

where W and K are the Newtonian potential and kinetic energy (per unit mass), 

W= l -{{l + S)cp), K=l{(l + 6)v 2 ). (42) 

Here, 5 = bpj (p) is the density perturbation in the matter rest frame and v is the peculiar velocity. The argument is 
that observationally both W and K are small, hence, inhomogeneities cannot mimic dark energy. 

However the model constructed in Sec. Mil seems to belie this simple and compelling argument. If one interprets 
the observations of p(z) of the LTB model in terms of the fiducial ACDM phenomenological background solution, one 
would conclude that the perturbations are very small (viz. zero), hence W and K zero, and inhomogeneities could 
not mimic dark energy. But, at least for di,(z) and p(z), inhomogeneities do exactly that! 
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We point out that the very notion of peculiar velocities requires a background solution (a norm) about which 
to define "peculiarity" [49j. Also, Enqvist, Mattsson, and Rigopoulos [5(| explicitly demonstrate how spherically 
symmetric perturbations of a dust Universe can lead to acceleration while the perturbations of the gravitational 
potential remain small on all scales. Perhaps the general lesson is that there are highly non-Gaussian inhomogeneities 
in the late Universe, and the coherence of structures causes small deviations in observables to sum to a large deviation. 

We can venture a different way to frame the argument: If the observed density fluctuations, 5p(z), originated from 
the growth of perturbations in an FLRW cosmology, and the nonlinear perturbations on one scale did not affect the 
growth of structure or the evolution of the background solution on any other scale, and inhomogeneities do not show 
coherence or non-Gaussianity, and <fi ~ ip, then one can say that the interpretation of observations would yield a value 
of W and K that are too small to mimic dark energy. 

This way of framing the argument makes manifest the fact that there are dynamical assumptions that enter, and it 
cannot simply be a matter that small peculiar velocities rule out the backreaction proposal. If any of the dynamical 
assumptions fail, the argument fails. 

Another lesson that can be extracted from the exercise of this paper concerns using averaging procedures to quantify 
the effect of inhomogeneities on observables. This was the subject of the previous section. The conclusion was that 
while the two averaging prescriptions discussed give an indication that the LTB model would have observables that 
might be interpreted as arising from dark energy, the effective equation of state parameters do not resemble that 
of the fiducial ACDM model. This might be because the averaging approach is not very useful, or it might imply 
something quite troubling. Perhaps the equation of state parameter iu, which has received so much attention, has 
a physical interpretation only in FLRW models, and in an inhomogeneous Universe, where there is no natural time 
slicing, it is an unphysical construct. 

This paper has developed and discussed issues related to the interpretation of cosmological observables in con- 
structing cosmological models and theories. Observations of di,(z) and p(z) do not, by and of themselves, prove that 
there is dark energy or that the Universe is accelerating in the usual sense; they do so only if one assumes that the 
Universe is homogeneous and isotropic and the dynamics of the Universe are governed by general relativity. This 
paper discussed aspects of the program to develop understanding of how an inhomogeneous Universe might mimic 
dark energy. While consideration of the results of this investigation might lead one to conclude that there are no 
arguments that conclusively exclude the backreaction proposal, we are no closer to developing an inhomogeneous 
cosmology alternative to dark energy. 

One attempt to construct inhomogeneous cosmological models involves studying "Swiss-cheese" models comprised of 
low density voids embedded in an Einstein-de Sitter background (HI [H [H, M, M, HE M, S [E| . If the void regions 
are packed together the regions between them form a network of high-density filamentary, or spaghetti-like, structures. 
Further refinements of this approach have included spherical mass concentrations referred to as "meatballs" (60| . A 
Swiss-cheese, meatball, and spaghetti cosmological model may include ingredients that are desirable on their own, or 
work well in combination with one other ingredient, but combining them all may lead to a model that is not easily 
digestible. 

Whatever direction the study of inhomogeneous cosmologies might take, consideration of the results of this inves- 
tigation should be important. 
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